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Within the approximation of dominant charging collisions 
the explicit microscopic calculations of the Fokker-Planck ki- 
netic coefficients for highly-charged grains moving in plasma 
are performed. It is shown that due to ion absorption by grain 
the friction coefficient can be negative and the appropriate 
threshold value of the grain charge is found. The stationary 
solutions of the Fokker-Planck equation with the velocity de- 
pendent kinetic coefficient are obtained and the considerable 
deviation of such solutions from the Maxwellian distribution 
is established. 

PACS number(s) : 52.27.Lw, 52.20.Hv, 52.25.Fi 

In the systems close to equilibrium Brownian parti- 
cles keep stationary random motion under the action of 
stochastic forces which are compensated by the particle 
friction and thus, the work produced by the Langevin 
sources is equal to the energy dissipated in course of 
the Brownian particle motion. This energy balance is 
described by the fluctuation-dissipation theorem in the 
form of the Einstein law. Obviously, the fluctuation- 
dissipation theorem and the Einstein relation can be vi- 
olated in the case of non-equilibrium systems (even in 
the stationary case), in particular in the open systems. 
Starting from the classical Lord Rayleigh work [1] many 
studies of the non-equilibrium motion of Brownian par- 
ticles with additional (inner, or external) energy supply 
have been performed. In particular, such studies are of 
great importance for physico-chemical [2, 3] and biolog- 
ical [4] systems in which non-equilibrium Brownian par- 
ticle motion is referred as the motion of active Brownian 
particles. Recently the dynamical and energetic aspects 
of motion for the active Brownian particles have been de- 
scribed on the basis of the Langevin equation and the ap- 
propriate Fokker-Planck equation [5, 6]. The possibility 
of negative friction (negative values of friction coefficient) 
of the Brownian particles was regarded, as a result of en- 
ergy pumping. For some phenomenological dependences 
of the friction coefficient as a function of the grain's ve- 
locity the one-particle stationary non-Maxwellian distri- 
bution function was found. 

The traditional formulations of the non-equilibrium 
Brownian motion are based on some phenomenological 
expressions for the friction and diffusion coefficients. In 
particular, it means that deviations from the Einstein 
relation, as well as the velocity dependence of these co- 
efficients are postulated and high level of uncertainty for 
the application of such models to the real systems takes 
place. Here we will consider another situation, when the 
kinetic coefficients can be calculated exactly on the ba- 
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sis of the microscopically derived Fokker-Planck equation 
for dusty plasmas [7,8] . It will be shown that in the case 
of strong interaction parameter T = e 2 Z g Zi/aTi ^> 1 
(here Z g , Zi are the charge numbers for the grains and 
ions respectively, a is the grain radius, Ti is the ion 
temperature) the negative friction coefficient appears for 
some velocity domain. If the charging collisions are 
dominant, this domain is determined by the inequal- 
ities: v 2 < io?4;(r 
(r - 1) « 1,(4 

The physical reason for manifestation of negative friction 
is clear: the cross-section for ion absorption by grain in- 
creases with the relative velocity between the ion and 
grain decrease due to the charge-dependent part of the 
cross-section. Therefore, for a moving highly-charged 
grain (r > 1) the momentum transfer from ions to the 
grain in the direction of grain velocity is higher than in 
the opposite direction. 

Naturally, the Coulomb scattering and particle friction 
related to the ion-grain and neutral-grain elastic colli- 
sions will increase the threshold for negative friction and 
even can suppress it for some plasma parameters. How- 
ever, the latter are rather sophisticated processes, and 
their theoretical description within various approxima- 
tions can shadow the physics of the phenomenon. Some 
estimate of the influence of such processes will be also 
done below, but the detailed description of their influ- 
ence should be a matter of further consideration. On 
the other hand, the problem of microscopic investigation 
of negative friction in dusty plasmas is so fundamental, 
that it deserves a description of its simplest manifestation 
(which occurs in the case of dominant charging collisions) 
to be considered in this Letter. 

We start from the Fokker-Planck equation for the 
spherical grains in dusty plasma with typical narrow 
charge distribution around a negative value q = e e Z g , 
that permits to put all grain charges to be equal. We 
also ignore the increase of the grain mass [9, 10] assum- 
ing that neutral atoms generated in the course of the 
surface electron-ion recombination escape from the grain 
surface into a plasma. Than, for the conditions ij 3> Ajj , 
a, where U is the ion mean free path length and is 
the plasma screening (Debye) length, the friction and dif- 
fusion coefficients in the Fokker-Planck kinetic equation 
/3(g,v) and D(q,v) [8] are given by: 

/3(?,v) = - —fdV^a ag ( q ,W-v>\)\v-V\ 
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Here cr SQ (<?, v) is the cross-section for grain charging 
within the orbital motion limited (OML) theory. In this 
approximation all electrons and ions approaching to the 
grain on the distance less than a are assumed to be ab- 
sorbed. Subscript a = e, i labels plasma particle species. 
The rest of notation is traditional. To include the pro- 
cesses of the electron, ion and atom scattering we have 
to summarize the appropriate coefficients on the different 
type of the processes. In the case of dominant charging 
collisions integration in Eq. (1) can be performed explic- 
itly. The ion part of f3 (which exceeds the electron one 
at least in (rriiT e / 'meTi) 1 / 2 times) is the following 
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Equation (4) has a root T](T), which exists and is small 
(according to the conditions of applicability for this ex- 
pansion) only for r > 1, but f-l<l: 



V (T) = 5 
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Therefore the function ^(77, T) is negative at 77 < 77(F), if 
r > 1. Eq. (4) shows that for T < 1, when /3j is positive 
for all 77, the derivative (d@i(j), r) / 0^77) | ?i= o changes its sign 
at r = 1/3. Near 77 = the friction (3i decreases as 
function 77, if T < 1/3 and increases, if T > 1/3. 




Here n a are the densities of the electrons and ions, Erfy 7 ?? 
is the error function and 77 = 77(77) = v 2 /2v 2 ?i . It is 
easy to see that the terms Ji(?7) and ^2(77, T) describe 
the parts of (3i(q,v) related to purely geometrical and 
charge-dependent collecting cross-sections, respectively. 
The integration for Di(q,v) leads to the expression: 
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As follows from the physical reason and directly from 
Eq. (3), the coefficient D(q,v) is always positive. At the 
same time for some values of 77 and T the friction coeffi- 
cient (3i(q,v) can be negative. To find the root 77(F) of 
the equation f3i(rj, T) — let us consider for two limiting 
cases 77 <C 1 and 77 ^> 1. For 77 <C 1 Eq. (2) gives 
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FIG. 1. Numerical solution of the equation /?(?/, T) = 
separating the positive and negative values of the friction co- 
efficient. 

For 77 3> 1 and arbitrary values of T Eq. (2) gives 
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It means that for T ^> 1 and large 77 the equation 
ft^T) = has a root 7/1 (r) ~ (l+2r)/2 - T. Naturally, 
for all r > 1 there exists appropriate 7y(r), which is the 
root of that equation. This conclusion is confirmed by the 
exact numerical calculations of 7?(r) (Fig. 1). Asymptot- 
ically (for 77 » max(l, T)) ^(77, V) tends to zero as \Zij~ 1 
and is positive for 77 > 77(F). For the case of the negative 
friction (r > 1) the maximum of the coefficient ^(77, T) 
is located at the point ?7 m (r) ^> 1: 
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For the diffusion coefficient Di(rj,T) the expansions for 
7/ <C 1 and 77 ^> 1 lead to 
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The typical behaviour of Pi(r], T) and Di(r], T) calculated 
numerically on the basis of Eqs. (2), (3) is shown in 
Figs. 2, 3. 
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FIG. 2. The velocity dependences of the dimensionless fric- 
tion coefficient f3(r),T) — Ii(r]) + 12(77, F) for the different val- 
ues of r. 




FIG. 3. The same dependences for the dimensionless diffu- 
sion coefficient D(q,v) = Ki{rj) + K2(rj,T). 

The stationary solution of the Fokkcr-Planck equation 
with the kinetic coefficients (2), (3) for the grain distri- 
bution function f g (q,v) is 
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where C is a constant, providing normalization 
j dvf g (q,v) = 1. The velocity dependence of this so- 
lution for different values of T is shown in Fig. 4. 

In order to get some analytical estimates let us con- 
sider the vicinity of the point T = 1. In such case the 
integration in Eq. (9) leads to the non-Maxwellian distri- 
bution function, which for T > 1 possesses a maximum 
at v ^ 0: 
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and D;/ is the cylindrical parabolic function. Eq. (10) is 
relevant, if the integral over 77 converges, what is valid 
for positive ~f(F > 0,472). At the same time, as was 
mentioned above, for the applicability of the expansions 
(4), (8) the inequality r < 1 + S with 5 <C 1 is required. 
It is clear from Eqs. (6), (8), that the asymptotic be- 
haviour of the distribution function is non-exponential 
/g( ? 7)l>)^oo ~ r]~( mg ' mi >, but for the values of T under 
consideration it is not essential for calculations of the av- 
erages, due to the rapid convergence of the integrals over 
77. In particular, the average kinetic energy K of grains 
is 
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that gives K = 2, 86(m g /m, i ) 1/2 7 n i for V = 1. Thus, the 
ion absorption by grains can lead to grain heating and 
grain average kinetic energy (or their effective temper- 
ature) can be much higher, than the electron and ion 
temperatures. This effect can be applied to the explana- 
tion of the experimental data [11, 12], as it was already 
suggested in [7,8] on the basis of the velocity independent 
approximation for p and D. 

Let us consider now the case r> 1. In this case the 
ratio Q = —(Ti/mi)0(r],T)/D(r],T) can be represented 
for all values of 77 with a good accuracy by the function 
Q = (m g /mi)(r - 77)/ [(1 + T)(?7 + r)j that gives 
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It is evident, that f g (r],T) given by Eq. (13) is non- 
exponential. 



3 



Finally for the domain T < 0, 472 we can omit the term 
~ v 4 in Eq. (10) and the distribution f g (q,v) becomes 
Maxwellian with the effective grain temperature: 
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The above conclusions are in a good agreement with 
the results of numerical calculations of / 9 (ry,r) on the 
basis of Eq. (9) (Fig. 4). 



In conclusion, ion absorption by grains can generate 
negative friction and provide the substantial increase of 
the average grain kinetic energy in comparison with the 
temperatures of the other plasma components. Micro- 
scopical justification of negative friction on kinetic level 
is presented and deviation of the grain distribution func- 
tion from the Maxwellian distribution is found. For dusty 
plasma, as for an open system, the fluctuation-dissipation 
theorem in the form of Einstein relation is not applicable. 
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FIG. 4. The same dependences for the distribution function 

Uv,r). 

If we take into account the processes of atom-grain 
and ion-grain elastic scattering the coefficients /3(r], T) 
and D(rj, T) should include additional terms. To esti- 
mate, for example, contribution of atom-grain and ion- 
grain scattering we can consider f3(rj) calculated on the 
basis of the Eq. (2) with the appropriate transport cross- 
sections. For the case rj <C 1 the friction coefficient is 
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The negative friction exists for small rj, if the Coulomb 
scattering is strongly suppressed, when the Coulomb log- 
arithm In A is small [13-15], that is typical for strong in- 
teraction. The root of the function P(rj, T) is shifted to 
the region of large, T, which is determined by the atom 
density n a . The result of rough estimate for the case 
Ali > A_d > a (Ali is the Landau length) gives 
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In general, a more detail consideration of scattering pro- 
cesses and the effects of strong interactions between the 
grains [15], is needed for the exact description of the re- 
gion of negative friction. 
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